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ABSTRACT 

Several interesting luminosity correlations among gamma-ray burst (GRB) 
variables have been recently discussed extensively. In this paper, we derive the 
six luminosity correlations (ri ag - L, V - L, E pc&k - L, E pc&k - E 7 , t rt - L, 
-E'peak — -E^.iso) from the light curves and spectra of the latest 116 long GRBs, 
including the time lag (r\ ag ) between low and high photon energy light curves, 
the variability (V) of the light curve, the peak energy of the spectrum (-E pe ak), 
and the minimum rise time (trt) of the peaks. We find that the intrinsic scat- 
ter of the V — L correlation is too large and there seems no inherent correlation 
between the two parameters using the latest GRB data. The other five correla- 
tions indeed exist when the sample is enlarged. The -E pca k — E 1 correlation has 
a significantly lower intrinsic scatter compared to the other correlations. We 
divide the full data into four redshift bins when testing possible evolution of 
the correlations with redshift. We find no statistically significant evidence for 
the redshift evolution of the luminosity correlations. To avoid the circularity 
problem when constraining the cosmological parameters, we simultaneously 
minimize x 2 with respect to both correlation parameters a, b and the cosmo- 
logical parameters using the maximum likelihood method. For the flat ACDM, 
the best fit is f2 m = 0.31^q \q. We also constrain the possible evolution of the 
equation of state (EOS) of the dark energy using the GRBs together with the 
Union2 compilation of SNe la and the H (z) data. The result is consistent with 
the cosmological constant at 2a confidence level and mainly due to the GRB 
data, the dark energy EOS shows slight deviation from -1 at z ^ 0.5 as was 
persistently presented with many previous data sets. 

Key words: cosmology: observations - gamma rays: bursts - cosmology: 
distance scale - cosmology: cosmological parameters 



1 INTRODUCTION 

Unexpected accelerating expansion of the universe was 
first discovered by observing type la supernovae (SNe 
la) (Riess et al. 1998; Perlmutter et al. 1999). Inde- 
pendent observations from baryonic acoustic oscillations 
(BAO) (Eisenstein et al. 2005; Percival et al. 2007), 
the anisotropy spectrum of cosmic microwave back- 
ground radiation (Komatsu et al. 2009) and the large 



* fayinwang@nju.edu.cn 
f qishill@gmail.com 
| dzg@nju.edu.cn 



scale structure data from large galaxy redshift surveys 
(Tegmark et al. 2006) have confirmed this surprising re- 
sult. This acceleration is commonly attributed to dark 
energy, which is the most mysterious problem in modern 
cosmology. Among parameters that describe the prop- 
erties of dark energy, the equation of state (EOS) is 
one of the most important. Whether and how it evolves 
with time is crucial in distinguishing different cosmo- 
logical models. A nearly model-independent approach 
in which uncorrelated estimates are made about discrete 
w(z) at different redshifts has been extensively discussed 
(Huterer & Cooray 2005; Riess et al. 2007; Sullivan et 
al. 2007; Qi, Wang & Lu 2008). 
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In order to measure the expansion history of our 
Universe, we need the Hubble diagram of standard can- 
dles. SNe la are the well known standard candles that 
have played an important role in constraining cosmo- 
logical parameters. Unfortunately, it is difficult to ob- 
serve SNe la at z > 1.7, even with excellent space based 
projects such as SNAP (Aldering et al. 2004) . They can- 
not provide any information on the cosmic expansion 
beyond redshift 1.7. With gamma-ray bursts (GRBs), 
we can access much higher redshifts. The high lumi- 
nosities of GRBs make them detectable out to the edge 
of the visible universe (Lamb & Reichart 2000; Ciardi 
& Loeb 2000; Bromm & Loeb 2002, 2006). The far- 
thest GRB observed hitherto is GRB 090423 at z = 8.2 
(Tanvir et al. 2009; Salvaterra et al. 2009). Schaefer 
(2007) complied 69 GRBs to make simultaneous uses 
of five luminosity relations, which are the correlations 
of Ti ag - L (Norris, Marani & Bonnell 2000), V — L 
(Fenimore & Ramirez-Ruiz 2000), -E pca k — L (Schaefer 
et al. 2003; Wei & Gao 2003), £ pca k - £ 7 (Ghirlanda et 
al. 2004a), and r RT - L (Schaefer 2007). Here the time 
lag (ti^) is the time shift between the hard and soft 
light curves, L is the peak luminosity of a GRB, the 
variability V of a burst denotes whether its light curve 
is spiky or smooth and it can be obtained by calculat- 
ing the normalized variance of an observed light curve 
around a smoothed version of that light curve (Fenimore 
& Ramirez- Ruiz 2000), £p C ak is the photon energy at 
which the vF v spectrum peaks, = (1 — cos 0j)E^ t i SO 
is the collimation-corrected energy of a GRB, and the 
minimum rise time (trt) in the gamma-ray light curve 
is the shortest time over which the light curve rises by 
half of the peak flux of the pulse. More recently, Yu et al. 
(2009) found that, for the three-dimensional (3D) lumi- 
nosity relations between the luminosity and an energy 
scale Speak and a timescale (ri ag or trt), the intrin- 
sic scatters are considerably smaller than those of cor- 
responding two-dimensional (2D) luminosity relations. 
Dainotti et al. (2008, 2010) and Qi & Lu (2010) found 
new correlations between the transition times of the X- 
ray light curve from exponential to power law and the 
X-ray luminosities at the transitions. After being cal- 
ibrated with luminosity relations, GRBs may be used 
as standard candles to provide information on cosmic 
expansion at high redshifts and, at the same time, to 
tighten the constraints on cosmic expansion at low red- 
shifts (Dai et al. 2004; Ghirlanda et al. 2004b; Friedman 
& Bloom 2005; Liang & Zhang 2005, 2006; Wang & Dai 
2006; Schaefer 2007; Wright 2007; Wang, Dai & Zhu 
2007; Wang 2008; Qi, Wang & Lu 2008a,b; Liang et 
al. 2008; Amati et al. 2008; Cardone et al. 2009, 2010; 
Liang et al. 2009; Qi, Lu & Wang 2009; Izzo et al. 
2009; Liang & Zhu 2010). GRBs also can potentially 
probe the cosmographic parameters to distinguish be- 
tween dark energy and modified gravity models (Wang, 
Dai & Qi 2009a, b; Vitagliano et al. 2010; Capozziello 
& Izzo 2008). 

The correlations among GRB variables span a very 
large range in redshift. Possible evolution effect must 
be considered when we use these correlations. Li (2007) 
used the Amati relation (-B pca k — E lt i so ) (Amati et 
al. 2002) as an example to test the cosmic evolution 



of GRBs and found that the slope of the correla- 
tion evolves with the redshift. In contrast, Basilakos 
& Perivolaropoulos (2008) found no statistically sig- 
nificant evidence for redshift dependence of correlation 
slopes using 69 GRBs. In this paper, we first enlarge 
the GRB sample with the new data from Xiao & Schae- 
fer (2009). Our sample includes 116 GRBs ranging from 
z = 0.17 to z = 8.2. We divide these GRBs into four 
redshift bins to investigate the possible evolution effect. 
Here the focus is on the correlations, so we fix the cosmo- 
logical parameters. We also use GRBs to constrain the 
cosmological parameters and dark energy EOS. In order 
to avoid the circularity problem, we simultaneously fit 
the correlation parameters and the cosmological param- 
eters. 

The structure of this paper is as follows: in the next 
section we show the latest GRB data and describe our 
fitting methods. In section 3 we present the updated lu- 
minosity correlations and test their redshift dependence. 
Constraints on cosmological parameters and equation of 
state of dark energy are presented in section 4. Some 
conclusions are presented in section 5. 



2 OBSERVATIONAL DATA AND 
ANALYSIS METHOD 

The luminosity correlations we will discuss here typi- 
cally relate a GRB observable with the isotropic peak 
luminosity L (it is also referenced to as L p in many 
papers), the isotropic energy -B 7l i so , or the collimation- 
corrected energy _E 7 . The isotropic peak luminosity is 
given by 



L = 47rd t P bo io, 

the isotropic energy is 

E~,a so = 47rd|S b oio(l + z) _1 , 

and the collimation-corrected energy is 



Ery E~> isO-^bc 



47Td t S'bolofb e am(l + A 



(1) 



(2) 



(3) 



Here, Pboio and Sbolo are t ne bolometric peak flux and 
fluence, respectively, while Pbeam = 1 — cos 9j et is the 
beaming factor. From Sari, Piran, & Halpern (1999), 

0jet = 0.161[W(1 + z )f S ( n ^/^7,iso,52) 1/S , (4) 

where z is the redshift, ij e t is the jet break time mea- 
sured in days, n is the density of the circumburst 
medium in particles per cubic centimeter, rj-y is the ra- 
diative efficiency, and -E 7 ,i so ,52 is the isotropic energy in 
units of 10 52 erg for an Earth- facing jet. The jet break 
time (tjet) can be measured when the afterglow bright- 
ness has a power-law decline that suddenly steepens due 
to the slowing down of the jet until the relativistic beam- 
ing roughly equals the jet opening angle. In the absence 
of these detailed fits, we adopt ry 7 = 0.2 and n = 3 cm -3 
(Schaefer 2007). Note that Pboio and Sbolo are computed 
from the observed GRB energy spectrum &(E) as fol- 
lows (Ghirlanda et al. 2004a, Amati 2006): 



P x 



*(E)dE 



(5) 
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5boio — 5* 



fi/(i+ 1 z ) Z ' > E<&{E)dE 
E$(E)dE 



(6) 



with P and S being the observed peak energy and flu- 
ence in units of photons/cm 2 /s and erg/cm 2 , respec- 
tively, and (-Emm, P max ) the detection thresholds of the 
observing instrument. For pre-Swift GRBs, we take the 
values of Pboio and SWo directly from Schaefer (2007). 
For those GRBs observed by Swift, we adopt the values 
of P and S from Swift website Q and calculate Pboio and 
Sboio using the above formulae. Concerning the errors of 
Pboio and Sboio during the calculation, we only take into 
account the errors propagating from that of P and S. 
The uncertainties from <&(E) are absorbed into intrinsic 
scatters of the correlations. Note that the energy spec- 
trum is modeled using a smoothly broken power -law 
(Band et al. 1993), 



$(P) = 



BE 13 



(7) 



otherwise 



where a is the asymptotic power-law index for photon 
energies below the break and /3 is the power-law index 
for photon energies above the break. We use the val- 
ues of a and ft from Xiao & Schaefer (2009). The lu- 
minosity correlations are power-law relations of either 
L, P 7 ,iso or P 7 as a function of Ti ag , V, P pca k, or t rt . 
The luminosity indicators of ri ag , V, P pca k, and trt 
are also directly taken from Xiao & Schaefer (2009). L, 
-Sy.iso, and E 1 depend not only on the GRB observables 
Pboio or 5boio, but also on the cosmological parameters 
through the luminosity distance cLl, which in a flat uni- 
verse is expressed in terms of the Hubble expansion rate 
H(z) = H E(z) as 



d L {Qm, z) 



(1 + ^Po 



dz' 



(8) 



where E 2 (z) = O m (l + zf + Q^f x (z) and the dmien- 
sionless dark energy density f x (z) is given by (w(z) is 
the EOS of dark energy) 



/ x (z) = exp 



1 + w(z) 
l + i 



dz 



(9) 



When the focus is on the luminosity correlations them- 
selves, the cosmological parameters here are fixed. 

The luminosity correlations involved in this paper 
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Concerning the luminosity indicators in the correla- 
tions, for the temporal indicators, the observed quan- 
tities must be divided by 1 + z to correct the time dila- 
tion. The observed V-value must be multiplied by 1 + z 
because it varies inversely with time, and the observed 
P pca k must be multiplied by 1 + z to correct the redshift 
dilation of the spectrum. 

The first five of the correlations listed above were 
the ones considered in Schaefer (2007) and Xiao & 
Schaefer (2009). We add in our analysis the investigation 
of the Amati correlation (P pea k — P 7 ,i S o), which was ini- 
tially discovered on a small sample of BeppoSAX GRBs 
with known redshift (Amati et al. 2002) and confirmed 
afterwards by Swift observations (Amati 2006). Com- 
pared to the Pp Ca k — E-y correlation, due to the indepen- 
dence of #j c t, the P poa k — p7,iso correlation can be used 
for almost the whole GRB sample and does not suffer 
from the assumptions and uncertainties around #j e t that 
affect the P pGa k — P 7 correlation. Also, compared to the 
Ep ea k — L correlation, the P pGa k — P 7 ,iso correlation is 
not affected by assumptions on the peak flux time scale 
and on the spectral shape at the peak (i.e., the peak 
luminosity is always computed by assuming the spec- 
tral shape of the time-averaged spectrum, which is not 
physical, given that the spectrum at the peak is often 
much different than the average one). There is also de- 
bate about the reality of these correlations, see Nakar & 
Piran (2005), Band & Preece 2005, Butler et al. (2007), 
Butler et al. (2009), Ghirlanda et al. (2005), Bosnjak et 
al. (2008), Ghirlanda et al. (2008), Nava et al. (2008), 
Krimm et al. (2009), Amati et al. (2009), Ghirlanda et 
al. (2010) etc. 

In Table 1, we list the variables of 116 GRBs that 
we use in fitting luminosity correlations. In addition to 
the GRBs included in the analysis of Xiao & Schaefer 
(2009), we add the GRB090423, which has the highest 
redshift so far. We use the typical spectral index a = — 1 
and P = —2.2 for this burst (Salvaterra et al. 2009). 



1 Sec http:/ /swift. gsfc.nasa.gov/docs/swift/archive/grb_tablc 
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61 
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± 


3 
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6 
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± 
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0.061 ±0.033 
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1 


01 ± 0.75 
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83 
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± 


1 


44E-08 
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-5 
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± 


1 
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1.12E-06 ± 1.19E-07 
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11 ± 2.28 
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3 


35 


9.83E-08 


± 


1 
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1.09E-06 ± 6.98E-08 
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41+? 
-5 


1 


10 ± 1.47 
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2 


20 


1.93E-06 


± 


5 


18E-08 


5.09E-06 ± 1.03E-07 










0.06 ± 0.01 


0.015 ±0.003 


198+?? 





13 ± 0.01 


051016B 





94 


1.92E-07 


± 


1 


43E-08 


4.31E-07 ± 3.39E-08 












0.008 ± 0.030 


24+! 
-7 






051022 





80 


1.1E-5 


± 


8 


7E-7 


3.40E-4 ± 1.2E-5 
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± 
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0.088 ± 0.008 
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2 


35 


8.30E-07 


± 


8 
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70 ± 1.25 
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1 
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± 
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3 
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1 


11 ± 1.71 
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4 
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4.41E-07 


± 


1 
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0.01 ±0.03 


0.007 ± 0.004 
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1 
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3 


91 


5.37E-07 


± 


3 


36E-08 


1.97E-05 ± 6.39E-07 










0.15 ±0.17 


0.183 ±0.033 


149+«° 





73 ± 0.50 


060223A 


4 


11 


2.06E-07 


± 


1 


67E-08 


1.51E-06 ± 6.54E-08 












0.036 ±0.021 


71+100 





41 ±0.23 


060418 


1 


49 


1.49E-06 


± 


1 


85E-08 


2.62E-05 ± 4.85E-07 










0.22 ±0.03 


0.104 ±0.008 


230+jO 





67 ±0.08 


060502A 


1 


51 


3.72E-07 


± 


2 


81E-08 


6.59E-06 ± 1.77E-07 










4.90 ±0.11 


0.004 ±0.010 


156+J°° 


2 


94 ± 1.19 


060510B 


1 


90 


9.51E-08 


± 


1 


12E-08 


9.98E-06 ± 2.62E-07 












0.110 ±0.060 


95if 
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F 
[erg 
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/cm' 
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[erg/ cm 2 ] 


-^bcam 




[sec] 


v 


060512 
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32E 
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08 


6.04E-07 


± 
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5 


11 
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± 
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45E 


08 


2.42E-06 


± 


1.43E-07 
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3 
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33E- 
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± 


1 


53E 
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± 
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±0.09 
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2 


68 
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08 


± 


1 
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± 
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3 


80 
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08 


± 


1 


36E 
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1.58E-06 


± 


1.24E-07 


8.23E-04 d 


z 5.14E-05 
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3 


08 


2 


66E 


07 


± 


1 
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1 
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3 
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07 


± 
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12E 
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± 


1.96E-07 
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2 


71 


2 


30E- 


07 


± 


1 


42E 


08 
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± 


2.47E-07 












0.125 ± 0.022 


060729 





54 


1 


93E- 


07 


± 


1 


30E 


08 


6.43E-06 


± 


3.16E-07 
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o 


84 


1 


83E 


06 


± 


1 


44E 


08 


4.94E-05 


± 


4.91E-07 









29 
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o 


70 
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07 


± 


2 
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± 


2.17E-07 









36 
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2 


13 


6 


69E 


07 


± 


3 


36E 
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± 
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26 
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060926 


3 


21 


1 


56E- 


07 


± 


1 


22E 


08 
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± 
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1 


03 
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0.148 ± 0.050 


060927 


5 


60 


1 


02E 


07 


± 


1 


54E 


08 


2.37E-06 


± 


8.67E-08 









12 


±0.04 


0.094 ± 0.010 


061007 


1 


26 


7 


20E- 


06 


± 


1 


HE- 


07 


2.24E-04 


± 


1.72E-06 









11 


±0.01 


0.066 ± 0.003 


061110A 





76 


9 


79E- 


08 


± 


1 


SSE 


08 


2.71E-06 


± 


1.18E-07 












-0.038 ± 0.050 


061110B 


3 


11 


1 


79E- 


07 


± 


2 


66E 


08 


6.12E-06 


± 


3.38E-07 









24 


±0.36 


0.155 ± 0.064 


061121 


1 


31 


8 


04E- 


06 


± 


1 


07E 


07 


6.53E-05 


± 


5.76E-07 









0?, 


±0.01 


0.050 ± 0.003 


061222B 


3 


36 


2 


29E- 


07 


± 


3 


15E 


08 


5.01E-06 


± 


2.49E-07 












0.024 ± 0.043 


070110 


2 


35 


1 


12E- 


07 


± 


1 


36E 


08 


4.04E-06 


± 


1.64E-07 












-0.010 ± 0.031 


070208 


1 


17 


1 


39E- 


07 


± 


2 


06E 


08 


1.06E-06 


± 


1.46E-07 












0.083 ± 0.211 


070318 





81 


4 
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07 


± 


2 
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08 


7.34E-06 


± 


2.01E-07 












0.037 ± 0.008 


070411 


2 


95 


1 
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07 


± 


1 


31E 


08 


6.29E-06 


± 
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070506 


2 


31 


1 


67E- 


07 


± 


1 


38E 


08 


5.16E-07 


± 
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2 


52 


±0.04 
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82 


7 


67E- 


06 
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1 


18E 
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± 


6.15E-07 
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± 


5 
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01 
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070529 


2 


50 


3 


32E- 


07 


± 


5 
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7.44E-06 


± 


4.31E-07 
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070611 


2 


01 
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45E- 


07 


± 


2 


25E 


08 


9.52E-07 


± 


8.44E-08 












0.053 ± 0.080 


070612A 





62 


2 


77E 


07 


± 


1 


24E 


08 
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± 
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070714B 





92 


3 


24E- 


06 


± 


1 


46E 
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± 


6.77E-07 









03 


±0.01 
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2 


45 


6 


38E- 


08 


± 


9 


69E 
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± 
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070810A 


2 
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77E 
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1 


77E 
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± 
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1 


09 
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1 


10 
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71E 


06 


± 


1 


82E 
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± 
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17E 
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95 


9 


20E- 
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18E 
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8.37E-06 
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8 


61E 
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83E 
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0.009 ± 0.003 
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0.038 ± 0.021 
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080319C 


1 


95 


2 


22E 


06 


± 


7 


79E 


08 


1.77E-05 


± 


2.99E-07 
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080413A 


2 


43 
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2 
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± 


1.71E-07 
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±0.03 
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1 


10 


3 


17E 


06 


± 


8 


25E 


08 


8.00E-06 


± 


1.52E-07 
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±0.01 
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15E 
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± 
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±0.08 
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3 


20 


2 


77E 


07 


± 


2 


80E- 
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5.88E-07 


± 


5.50E-08 









15 
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0.168 ±0.055 
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1 


55 


8 


23E- 
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0.037 ±0.098 


080603B 


2 


69 


7 


57E 


07 


± 


2 


63E- 


08 


7.02E-06 


± 


1.78E-07 









08 
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3 


04 


8 
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2 
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080707 


1 


23 


1 


68E- 


07 


± 


1 


02E 


08 
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8.87E-08 
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080721 


2 


60 


9 


57E 


06 


± 


5 
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07 
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13 


±0.05 
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090423 
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2 


17E- 


07 


± 


1 
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-^pcak 

[keV] 



22r 6 

80ij2 2 
25^ 

40 ±5 
901S 
120±™° 
63l« 
103±?J 

61±1 

80t™ 

2011J 
72l 7 5 
399±" 
90±H 

M 7+ ; >3 
0i '-53 

606±^ 

49±g 
110 +30 

83lH 

311? 
233ly 
22211^ 

I8OI52 
92130 

871- 
1120+ 473 

,u -25 
44+9 

799lgj 
97+10 

52lg 
24± 7 7 
278lif 

7 o + 30 
'°-30 
7 o + 15 

28l« 
65ll 8 8 
3071^ 

20l 9 9 

259l 2 i 

170±« 
70+10 
'■^-io 

80±ig 
66l 2 * 

12+4 

24611? 

394+H 
70+20 

'°-20 

485l« 

48.6l^; 8 g 



TRT 

[sec] 



0.38 ±0.11 

1.22 ±0.72 

1.23 ±0.68 



1 


65 


± 


0.24 


1 


00 


± 


0.16 





52 


± 


0.09 



0.46 ±0.12 
0.38 ±0.02 

0.79 ±0.64 
0.98 ±0.19 



0.72 ±0.24 






12 


± 


0.06 





20 


± 


0.01 





58 


± 


0.06 



2.49 ± 1.48 
0.45 ± 0.04 

0.73 ±0.22 
0.88 ±0.07 

1.21 ±0.03 

0.20 ± 0.02 






57 


± 


0.44 





11 


± 


0.55 





14 


± 


0.01 





21 


± 


0.12 






65 


± 


0.01 





23 


± 


0.03 





50 


± 


0.03 





76 


± 


0.12 






22 


± 


0.03 





22 


± 


0.01 





18 


± 


0.06 



0.09 ± 0.04 



Table 1: The data of 116 GRBs used in our analysis. For pve-Swift GRBs, we 
take the values of Pbolo an d Sboio directly from Schaefer (2007). For those 
GRBs observed by Swift, we adopt the values of P and S from Swift website 
and calculate Pbolo an d Sbolo- We use the -Fboam value from Ghirlanda et al. 
(2007). Other data are taken from Xiao & Schaefer (2009). 
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The six luminosity relations can be expressed, in 
general, as R = AQ b and Eq. (|10[) - (|15[ ) are the corre- 
sponding logarithm forms 



log R = log A + b log Q =>• y = a + bx. 



(16) 



For the fit of this linear relation, we used the techniques 
presented in D'Agostini (2005), according to which, the 
joint likelihood function for the coefficients a and b and 
the intrinsic scatter 0" hit is 



L(a, b, (Ti n t) <x Y 



1 



x exp[ 



(Vi - 



+ < + b 2 c 



bxif 



2«t + < + 



(17) 



where Xi and yi are corresponding observational data 
for the ith GRB. When considering error propagation 
from a quantity, say £ with error ov, to its logarithm, 



we set 



iog(e+o-^)+iog(€-<T e ) 



and 



log(e + <T+)-log(f-<T f ) 



2 2 

the center value and the error of the logarithm corre- 
spondingly. This requires £ > er^~ (the quantities we are 
interested in here are all positive). Due to the limitation 
of the data, for a given luminosity correlation, not all 
the GRBs have all of the needed observational quanti- 
ties available and satisfy £ > a^ at the same time. The 
numbers of GRBs for each fit of the luminosity correla- 
tions are included in Table [2] 



3 TEST OF THE UPDATED LUMINOSITY 
CORRELATIONS 

3.1 Luminosity correlations 

Our fitting results for the six luminosity correlations 
are shown in Figure [1] and the last column of Table [2] 
We assume a flat ACDM with f) m = 0.27 and H = 
70 km s _1 Mpc _1 obtained from the five years WMAP 
data (Komatsu et al. 2009). The best-fit line and 2a 
confidence region are plotted in Figure [T] 

From Figure [T] we can see that the -E poa k — E-, 
correlation is the tightest one. The V — L relation is 
quite scattered. Its intrinsic scatter (a- m t = 0.67) has 
been larger than the one that could be expected for a 
linear relation. 



3.2 Test redshift variation of correlations 

In order to test if the correlations discussed in the above 
section vary with redshift, we divide the GRB samples 
into four groups corresponding to the following redshift 
bins: z 6 [0, 1], z 6 (1,2], 2 6 (2,3] and 2 G (3,8.5]. For 
each correlation and each redshift bin, we perform the 
same fit procedure as applied to the whole GRB sample 
to determine the parameters a, b and the intrinsic scat- 
ter <7int • The results of the fits and the number of GRBs 
used in each fit are summarized in Table [2] 

For further analysis, we perform linear fits to the 
parameters a and b versus redshift (the redshifts for the 
parameters a and b are calculated just by averaging the 
redshifts of the GRBs used in deriving corresponding a 
and b). These fits are shown in Figure [2] and the slopes 



of a and b versus redshift are presented in Table [3] For 
the -Epcak — E 1 correlation, there are no enough GRB 
samples to perform such fits. For the other luminosity 
correlations except for i? pca k — L, the slopes of b versus 
redshift are all consistent with zero at the 2a confidence 
level, and even for i5 pca k — L correlation, zero is near the 
edge of the 2a confidence interval of the slope of b versus 
redshift. Considering that, for the redshifts correspond- 
ing to a and 6, we only loosely use the average values of 
the redshifts of corresponding GRBs and the uncertain- 
ties in the redshifts are not taken into account (which 
leads to an underestimate of the uncertainties in the 
slopes), we can conclude that there is no statistically 
significant evidence for the evolution of the luminosity 
correlations with redshift. We didn't take into account 
the redshift evolution of the parameters a when drawing 
the conclusion, since they correspond to normalization 
factors in the luminosity correlations and a small change 
in b may lead to a larger change in a. In fact, as can be 
seen from Table [3l the slopes of a versus redshift for the 
first three luminosity correlations considerably deviate 
from zero. 



4 CONSTRAINTS ON COSMOLOGICAL 
PARAMETERS AND DARK ENERGY 
EQUATION OF STATE 

To constrain the cosmological parameters, we simulta- 
neously fit correlation parameters of GRBs and cos- 
mological parameters to avoid the circularity problem. 
Since the luminosity correlations of -E pca k — E~, and 
£pc a k — E lt i so describe almost the same physics, we can 
only include one of them in the fit to avoid strong cor- 
relation among the luminosity correlations. We choose 
the -E pca k — Ery correlation, which has a smaller intrinsic 
scatter. Since the intrinsic scatter of the V — L corre- 
lation has been too large, including it in the fit or not 
has little effect on the result. For the flat ACDM model, 
the combination of the correlations gives the result of 
fi m = 0.31 jl^'io ■ Schaefer (2007) used a combination 
of the same correlations with a smaller sample of GRBs 
and got the result of S7 m ~ 0.39. Our result is consistent 
with the value of Schaefer (2007) at the la confidence 
level. 

We also constrain the dark energy EOS using the 
GRBs together with SNe la and the H(z) data. We 
adopt the redshift binned parametrization for the dark 
energy EOS, as proposed in Huterer & Cooray (2005), in 
which the redshifts are divided into several bins and the 
dark energy EOS is taken to be constant in each redshift 
bin but can vary from bin to bin. For this parametriza- 
tion, f(z) = pde(z) I Pde(0) takes the form (Sullivan et 
al. 2007) 



< Z < Z n ) = (l+z) 3 ' 1 ^"' T[(l+ Zi ) 



i+l) 



(18) 



where Wi is the EOS parameter in the i th redshift bin 
defined by an upper boundary at Zi, and the zeroth bin 
is defined as zo = 0. Such a parametrization scheme 
assumes less about the nature of the dark energy, es- 
pecially at high redshift, compared with other simple 
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parametrizations, since independent parameters are in- 
troduced in every redshift range and it could, in princi- 
ple, approach any functional form with the increase of 
the number of redshift bins (of course, we would need 
enough observational data to constrain all the param- 
eters well). For a given set of observational data, the 
parameters Wi are usually correlated with each other, 
i.e. the covariance matrix 

C = (ww T ) - (w)(w T ), (19) 

is not diagonal. A new set of dark energy EOS param- 
eters wl defined by 

w = Tw. (20) 



is introduced to diagonalize the covariance matrix. The 
transformation of T advocated by Huterer & Cooray 
(2005) has the advantage that the weights (rows of T) 
are positive almost everywhere and localized in redshift 
fairly well, so the uncorrelated EOS parameters wl axe 
easy to interpret intuitively. The evolution of the dark 
energy with respect to the redshift could be estimated 
from these decorrelated EOS parameters. The transfor- 
mation of T is determined as follows. First, we define 
the Fisher matrix 

F ee C 1 = O t AO, (21) 

and then the transformation matrix T is given by 

T = O t A5Q, (22) 
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Correlation 



6 [0,1] 



e[i,2] 



z € [2, 3] 



z € [3,8.5] 



Total 



Tint 

GRB Number 



= 51.78 ±0.13 
= -0.55 ±0.17 
0.47 ±0.11 
15 



= 52.47 ±0.08 
= -0.77 ±0.13 
0.34 ±0.07 
19 



a = 52.44 ±0.12 
b= -0.83 ±0.21 
0.36 ±0.12 
11 



a = 52.57 ±0.24 
b= -0.60 ±0.38 
0.56 ±0.26 
8 



a = 52.28 ±0.07 
b= -0.77 ±0.10 
0.48 ± 0.05 
53 



V -L 



GRB Number 



a = 51.53 ±0.27 
b = 0.50 ± 0.44 
0.78 ±0.13 
26 



a = 52.25 ±0.24 
b = 0.35 ±0.26 
0.62 ±0.10 
25 



a = 52.38 ±0.48 
b = 0.25 ±0.40 
0.64 ±0.15 
15 



a = 53.03 ±0.42 
b= -0.23 ±0.33 
0.50 ±0.11 
15 



a = 51.86 ±0.15 
b = 0.59 ±0.15 
0.68 ±0.06 
81 



Speak — L 



GRB Numbcr 



a = 51.75 ±0.12 
b = 1.35 ±0.23 
0.63 ±0.09 
34 



a = 52.21 ±0.10 
b = 1.29 ±0.21 
0.48 ± 0.08 
30 



a = 52.27 ±0.09 
b = 1.40 ±0.24 
0.40 ± 0.07 
27 



a = 52.49 ±0.10 
b = 0.55 ±0.26 
0.43 ± 0.07 
25 



= 52.11 ±0.05 
= 1.40 ±0.12 
0.54 ±0.04 
116 



^pcak -C*7 



GRB Number 



a = 50.59 ±0.10 
b = 1.54 ±0.21 
0.24 ±0.11 
10 



= 50.66 ±0.11 
= 1.55 ±0.34 
0.14 ±0.13 
7 



= 50.60 ±0.07 
= 1.47 ±0.20 
0.31 ±0.08 
24 



TRT — L 

""int 
GRB Number 



= 52.49 ±0.18 
= -1.20 ±0.26 
0.50 ±0.08 
25 



a = 52.81 ±0.11 
b= -0.77 ±0.21 
0.39 ±0.07 
22 



a = 52.56 ±0.18 
b= -1.03 ±0.54 
0.55 ±0.16 
13 



a = 52.86 ±0.16 
b= -0.96 ±0.57 
0.43 ±0.14 
12 



a = 52.68 ±0.07 
b= -1.12 ±0.14 
0.48 ± 0.05 
72 



-^pcak -B"/ ,iso 



GRB Number 



a = 52.48 ±0.13 
b = 1.47 ±0.25 
0.62 ±0.10 
28 



a = 52.84 ±0.11 
b = 1.54 ±0.26 
0.54 ±0.09 
26 



a = 52.74 ±0.08 
b = 1.29 ±0.21 
0.36 ±0.07 
25 



a = 52.92 ±0.10 
b = 0.89 ±0.25 
0.39 ±0.07 
22 



a = 52.71 ±0.05 
b = 1.47 ±0.12 
0.49 ± 0.04 
101 



Table 2. Results of fits to the luminosity correlations for GRBs in each redshift bin and the whole sample. 





T lag — L 


V — L 


^poak — L Speak — E-y 


TRT — L 


-^pcak -^7,iso 


da/dz 


0.22 ±0.08 


0.44 ±0.15 


0.18 ±0.04 


0.06 ±0.07 


0.09 ± 0.04 


db/dz 


-0.06 ±0.12 


-0.23 ±0.15 


-0.22 ±0.10 


0.07 ±0.18 


-0.19 ± 0.10 



Table 3. The slopes of the parameters a and b versus redshift. 



except that the rows of the matrix T are normalized 
such that 



(23) 



We divided redshifts at points z = 0.2, 0.5, 1 and 
Markov chain Monte Carlo techniques are used with 
O(10 6 ) samples generated for each result. Since current 
observational data have only very weak constraints on 
the nature of dark energy at redshifts z > 1 (we tried 
constraining the dark energy EOS without imposing any 
prior on w(z > 1) using the parameterization described 
above with GRBs and other data sets, no substantial 
constraints on the dark energy EOS at redshifts z > 1 
can be obtained), we simply set w(z > 1) = —1, and 
focus on the dark energy EOS at z ^ 1. 

In addition to GRBs, we have used Union2 com- 
pilation of SNe la from Amanullah et al. (2010), BAO 



measurement from Percival et al. (2010) and Q. m h = 
0.213 ± 0.023 from Tegmark et al. (2004). We assumed 
the prior Q k = -0.014 ± 0.017 (Spergel et al. 2007) for 
the cosmic curvature. We also used the H(z) data from 
Stern et al. (2010) and Riess et al. (2009). 

For each luminosity correlation for GRBs, the Xgrb 
is calculated by 



2 

Xgrb 



-21nL, 



(24) 



where L is given by Eq. ()17p except that cosmological 
parameters are free parameters now. For other data set 
as well as the priors, the usual way of calculating x 2 is 
used, i.e., for a physical quantity £ with experimentally 
measured value £ , standard deviation <j£, and theoreti- 
cally predicted value £t(8), where 6 is a collection of pa- 
rameters needed to calculate the theoretical value, the 
X 2 value is given by 
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Figure 2. The correlation coefficients a, b obtained in four redshift bins for the five correlations. The -E pea k — ^7 correlation 
are not included because there are no enough GRBs in each redshift bin for this correlation. 



XW) = ■ (25) 

a c 

The total Xtotai is the sum of all the x 2 s from indepen- 
dent data. 

Figure [3] shows the result derived from the data 
set described above. We can see that though the dark 
energy is consistent with the cosmological constant 
(w(z) = —1) at the 2a confidence level, there is still 



considerable room for an evolving dark energy EOS. 
Notably, the slight deviation of the dark energy from 
the cosmological constant at z ^ 0.5, which persistently 
appears with many previous data sets, still exists here. 
For our result here, the deviation is mainly due to the 
GRBs. Though the Union compilation of SNe la gives 
the same trend of deviation from the cosmological con- 
stant due to the unexpected brightness of the Hubble 
data at z > 1, when the sample is enlarged, it seems 
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1.5 
1 

0.5 


13 -0.5 
-1 
-1.5 

-2 



-2,5 i ' ' ' ' ' 1 

0.2 0.4 0.6 0.8 1 1.2 

z 




1 . . 

0.8 - 

0.C - 

I 0.4 - 

0.2 

I I 

-0.2 1 

0.2 0.4 0.6 0.8 1 1.2 

Figure 3. Estimates of the uncorrelated dark energy EOS 
parameters ffij. Top panel: uncorrelated dark energy param- 
eters versus redshift, in which the vertical errorbars corre- 
spond to lcr and 2a confidence levels of and the horizon- 
tal errorbars span the corresponding redshift bins from which 
the contributions to Wi come most. Middle panel: Probabil- 
ity distribution for Wi. Bottom panel: window functions for 
Wi. 



that such a character of the SN la data has been aver- 
aged out in the Union2 compilation. As a comparison, 
we also present in Figure [4] the result derived from the 
data set without GRBs included. See also Wang et al. 
(2010) and Park et al. (2010) for similar analysis on the 
nature of the dark energy with Union2 compilation of 
SNe la. The deviation of the dark energy EOS from —1 
may arise from many possible reasons, for example, the 
statistical errors due to the limitation of current obser- 



1 

0.5 


-05 

-2.5 



0.2 0.4 0.6 0.8 1 1.2 

z 




1 

0.8 
0.6 



0.2 : 

I I 

-0.2 - 

0.2 0.4 0.C 0.8 1 1.2 

Figure 4. Estimates of the uncorrelated dark energy EOS 
parameters Wi. Same as the Figure [3] except that GRBs are 
not included in the fit. 

vational data, some biasing systematic errors in the ob- 
servational data (especially there is still some distance 
to calibrating GRBs as ideal standard candles), or the 
nature of the dark energy itself, etc. It should be made 
clear about the cause of the deviation with future ob- 
servational data in order to understand the dark energy 
better. 



5 CONCLUSIONS 

To build up the Hubble diagram to a redshift higher 
than the one of SNe la, most attempts have been de- 
voted to search for a method to make GRBs standardiz- 
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able candles. Different correlations have been proposed 
in order to build up a GRB Hubble diagram and con- 
strain cosmological parameters. As a further step, we 
have here considered the latest GRB dataset and lumi- 
nosity correlations to constrain the cosmological param- 
eters and dark energy. 

In this paper, we derived the six luminosity corre- 
lations (Ti ag — L,V — L, -Ep C ak - L, £ poak - E^, t rt - L, 
-Speak — -E 7l i S o) from the light curves and spectra of the 
latest 116 long GRBs. We find that the intrinsic scat- 
ter of V — L correlation is too large and there seems no 
inherent correlation between the two parameters using 
the latest GRB data. The other five correlations indeed 
exist when enlarging the sample. We have found no sta- 
tistically significant evidence for the redshift evolution 
of the luminosity correlations. However, even the best 
GRB luminosity correlation is currently not competi- 
tive with other cosmological probes of the cosmic accel- 
eration expansion since the cosmological parameter lcr 
errors derived from GRBs (fi m = CUl^io) 

are more 

than an order of magnitude larger than the correspond- 
ing errors obtained using SN la standard candles and 
other geometrical probes. But the estimates of cosmo- 
logical parameters from GRBs are important because 
they provide an independent confirmation of the results 
from other probes. 

We also performed an investigation on the dark en- 
ergy EOS using the GRBs together with the Union2 
compilation of SNe la and the H(z) data. The result 
is consistent with the cosmological constant at 2a con- 
fidence level. However, mainly due to the GRB data, 
the slight deviation of the dark energy EOS from —1 at 
2 ^ 0.5, which persistently appears with many previous 
data sets, still exists. 
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